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Abstract 

Basing on t — J model we calculate the k-dependence of a single hole photo- 
production probability for Cu02 plane at zero doping. We also discuss the 
radiation of spin-waves which can substantially deform the shape of photoe- 
mission spectra. 
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Recent photoemission measurement s3 for insulating Copper Oxide S^CuC^C^ give an 
unique possibility to determine single hole dispersion from experiment. However, to compare 
with data one needs to calculate the probability of the hole creation at a given momentum 
and energy transfer. The purpose of the present work is to calculate this probability in the 
framework of t — J model and elucidate the related question with momentum dependence 
of quasihole residue (see Refs.ii). In the conclusion we also comment on the shape of 
photoemission spectra. 

It is widely accepted that low energy dynamics of Cu0 2 planes in Copper Oxides is 
described by t — J model (see Ref.i for a review). At zero doping which corresponds to half 
filling this model is equivalent to Heisenberg antiferromagnet with long range Neel order 
in the ground state. We are interested in the process of an external photon kicking out 
an electron from the plane and creating a hole. The properties of a single hole are well 
established. For its description at t/J < 5 one can use selfconsistent Born approximation!. 
This approximation is justified by the absence of a single loop correction to the spin- wave 
vertexihi. In some sence it is analog of the well known Migdal theorem for electron-phonon 
interaction. Hole dispersion has minima at k = (±7r/2, ±7r/2) and it is almost degenerate 
along the face of magnetic Brillouin zone 7k = ^(cos k x + cos k y ) > 0. The hole wave function 
^kcr is characterized by quasimomentum k defined inside the magnetic Brillouin zone, and 
pseudo-spin a =][. The pseudo-spin denotes the sub lattice at which the hole is centered, 
and it is different from the usual spin. It is convenient to write the hole wave function in 
the form 

Vk. = 4Jo>, (i) 

where |0) is the ground state of the Heisenberg model, and /i kcr is the creation operator of 
the composite hole. Let us also denote by d ns the annihilation operator of an electron at 
the cite n and with usual spin s =ti- Let us fix the pseudo-spin a in Eq. (|l|), for example 
a =f. The quasiparticle residue of the composite hole is given by 
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We stress that the summation in this equation is restricted by one sublattice because the 
quasispin is fixed. This is the residue defined in the Refs. and due to the Bloch theorem 
the exact relation 
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is valid for translation at the inverse vector of the magnetic sublattice, (Q = (±7r, ±7r)). 

When a photon kicks out an electron from the system it does not separate the sublattices, 
therfore the external perturbation is 
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and the probability of the process is proportional to 
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Here N is number of cites in the lattice. Normalization of the external perturbation is rather 
arbitrary We choose it in such a way that v4 k = 1 for Ising case and t = 0. The value defined 
by Eq.(|5|) is different from the residue (0). For the calculation of v4 k it is convenient to use 
I = 1 string variational ansatz for the hole wave function!^ 
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where J2s denotes the summation over nearest neighbours, and the coefficients are given by 
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S k = (A 2 /4 + 4t 2 [(l + y )-(x + y) 7k ]) 
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We set J —I. The parameters A = 1.33, x = 0.56, y = 0.14, X = 0.8, Y = 0.72, u = 0.42, 
v — 0.12 are some combinations of the Heisenberg model spin-spin correlators calculated in 
the Ref.0. The residue (|2]) is proportional to the coefficient z/£ in the wave function (|6|). We 
would like to stress that the string representation of the hole wave function is dual to the 
usual spin-wave representation, but the physical matrix elements are certainly independent 
of the representation (see below). 

Keeping in mind the relations d n ^d n ^ = | + S z and dLd n f = S~ one can easily calculate 
the probability of photoemission (|^) in the minimal string approximation 
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Here a = | {0 1 jS^ 1 0) | = 0.3 is average magnetization, and qi = ~(0\S+ S~ +s \0) = —0.08 is the 
nearest neighbour spin-spin correlator for the Heisenberg model. Due to the 7k term in (H) 
^k+Q 7^ ^-k ' At zero hopping (t = 0) = ~ + a = 0.8. It is smaller than 1 because 
of spin quantum fluctuations in the initial state. This value agrees very well with result of 
spin-wave calculation^ and finite lattice computational. 

The / = 1 string ansatz (|6|) is a good approximation for the hole wave function at 
t < 1. For t > 1 further dressing by / > 1 strings, or, in other words, multy spin- wave 
virtual excitations is important. There are many multy spin- wave components in the wave 
function of a dressed hole. The assumption that these components do not give a coherent 
contribution in A k , but only reduce the probability of /i k j configuration is quite reasonable. 
This is exactly what happens for the hole spin-wave vertex§. Therefore 

A k = Z k 4 1} , (9) 

where Z k is the weight of in the exact /i kcr . The values of Z k =(7r/2,7r/2) ~ Z lk ~ calculated 
in the Ref.i for different t are presented in the Table I. Note that the quasiparticle residue 
(§) equals to Z k = Z k z/ k (l/2 + a),. Using this equality one can also find Z k from the results 
of Refs.iHi where quasiparticle residue 2 k was calculated. In Table I we present the ratio 
Aq/A = Aq /A^ for different values of t. We see that k-dependence qualitatively agrees 
with experimental datJ. We assume that t/J>0. It is interesting that at negative t the 
ratio Aq/Aq is inverted. At t < 1 our result for Aq/Aq reasonably agrees with that found 
by finite lattice computationslll'll! and numerical spin-wave analysis!. However at t — 3 the 
ratio Aq/Aq found in Refs.Bli! for finite lattices is bigger than that from Table I. Possible 
reason of this this disagreement is a variational nature of / = 1 string calculation. Therefore 
below we discuss more regular way of calculation of A^. 

Let us use not string, but spin-wave picture considering hopping term by pertur- 
bation theory! - !. The wave function of a bare hole (i.e. at t — 0) is of the form 
V'kT ~ \f^J2ne^ d n -^e l]irn \0). External perturbation is given by Eq. (|), and we can in- 
troduce bare vertex of single hole production. 

v£ = (^\M £4^10)^1. (io) 
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Here we neglect the spin quantum fluctuations in the initial state. External perturbation 
(|j) can also produce hole + spin-wave final state. Let us denote by the creation operator 
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of the spin- wave with momentum q and projection of spin S z = —1 (see Ref.0 for review) 
The vertex of hole + spin-wave creation equals 
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4(°l« q E^-|o) = W^ q; (ii) 



where t> q is Bogoliubov parameter dioganalizing spin-wave Hamiltonian: w q = J-^- + |, 
t> q = —sign^qjJ-^ — |, c<j q = 2^1 — t|, (see Ref.0) We stress that ( pTT| ) is a bare vertex. 
It corresponds to instant production of hole + spin wave, but not production of hole with 
subsiquent decay into hole + spin-wave. Note that Vfi sw -»• oo at q -»• 0. The reason is that 
perturbation (||) does not correspond any quasiparticle of the system, and therfore usual 
Goldstone theorem is not applicable. The verices ( |10D and (|TT| ) are presented at Fig. 1. 



Cross corresponds to external perturbation, solid line - to the hole, and dashed line - to the 
spin- wave. 

Bare hole-spin wave vertex is equal to g Piq = y^2/ (7 P M q + 7k u q ), with f = It (see 
e.g. Refs.iHi). It corresponds to the decay of the hole with momentum k into a hole with 
momentum p = k — q and a spin-wave with momentum q. 

Now we can easily calculate the first correction to the production vertex (|T0|). It is given 
by the diagram presented at Fig. 2. 
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This gives 



A k = (V b h + 6Vl 1] f = (1 + 0.4 f 7k ) 2 . (13) 



It is in a very good agreement with string variational result (H) which gives = 0.8(1 + 
0.37 t 7 k ) 2 at t -»• 0. 

To perform spin-wave calculations at large t one has to remember that one loop cor- 
rection to the hole-spin-wave vertex is absent§~! and therefore selfconsistent Born (SB) 
approximation! is valid. It means that only modification in comparison with naive pertur- 
bation theory is that in the diagram Fig. 2 one has to use dressed hole Green function found 
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in selfconsistent Born approximation. We carried out simplified computation of this type 
with dressed hole Green function in the form Z/(e — + zO). This is equivalent to the re- 
placement of bare hole-spin-wave coupling constant f = 2t with the effective one / = Z ■ f , 
see Ref.i. The results of this computation at t < 4 reasonably agree with variational answer 
®5 ©• One can certainly calculate more accurately using hole Green function found 
in Refs.i~i. However, we would like to stress that for detailed analysis of photoemission 
spectra is not enough. 

As we already discussed an external photon can produce hole + spin-wave as well as a sin- 
gle hole. Therefore the e, k dependence of photoemission spectrum is given by the imaginary 
part of diagrams presented at Fig. 3, with solid line be hole Green function in selfconsistent 
Born approximation Gsb- The Fig. 3 represents new Green function corresponding to the 
external perturbation (|J). The diagram Fig.3a gives usual Im Gsb(c, k). The diagram Fig. 
3b takes into account already discussed difference of A^ from pure residue Z^, but still this 
contribution is proportional to Im Gsb- The diagram Fig. 3c gives nonresonant background 
caused by the emission of spin-waves. And finally the diagrams Fig3d,e represent inter- 
ference between background and resonant production. This interference shifts maximum 
of photoemissin spectrum from e = ek- Whether this shift could explain the experimental 
speetrj is still an open question. 

Before conclusion we would like to note that the "electron" creation operator is actually 
Zhang and Rice singlet creation operator.0. Therfore formfactor of this singlet should 
be also taken into account. Calculation of F^ is not a problem because the wave function of 
the singlet is now well known@. 

In the present work we calculate single hole photoproduction probability A^ in frame- 
works of t — J model. In agreement with experiment A^ drops down outside magnetic 
Brillouin zone. We also point out that shape of photoemission spectrum is different from 
Im Gsb(c, k), and due to the interference the difference is not just a smooth background. 

I am very grateful to D. Poilblanc and G. Sawatzky for attracting my attention to this 
problem. This work has been done during my stay at the Laboratoire de Physique Quantique, 
Universite Paul Sabatier. I am gratefully acknowledge the hospitality and financial support. 
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TABLES 



TABLE I. The weight Z 7k «o of I = 1 string configuration in the dressed hole wave function, 
and the ratio of a single hole creation probabilities at k = Q = (±7r, ±7r) and k = for different 
values of t/J 



t/J 





0.2 


0.5 
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3 
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0.96 


0.88 


0.79 


0.63 


0.53 


0.45 


Aq/A 


1 


1.3 


1.9 


2.5 


3.3 


3.7 


3.9 



Figure captions, 

Fig.l Bare vertices of a single hole creation and hole + spin- wave creation. The cross 
corresponds to external perturbation, solid line - to the hole, and dashed line - to the spin- 
wave. 

Fig. 2 First order correction to the single hole creation amplitude. 

Fig. 3 The diagrams contributing into photoemission spectrum. Solid line corresponds to the 
hole Green function in selfconsistent Born approximation.. 
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